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I. INTRODUCTION 

The CERN Large Hadron Collider (LHC) provides the opportunity to study experimentally the production of scalar, 
pseudo-scalar, and vector mesons in peripheral collision of heavy ions. Peripheral kinematics implies the detection of 
\Q | particles produced in directions close to the axis of the colliding beams [l| . The main feature of these processes is the 
C^) . non-decreasing total and differential cross sections. The invariant mass of the created particles is assumed to be small 
in the fragmentation and central regions in comparison with the total energy in the center of mass of the colliding 
beams = 2£\ The application of the known theoretical approaches, such as the Nambu-Iona-Lasinio model as well 
as chiral perturbation theory, seems to be legitimate. We shall consider, in pure quantum electrodynamics (QED), the 
processes of the production of a single vector meson and of two vector mesons separated by a rapidity gap. It is known 
that the main contributions to the amplitudes of peripheral processes arise from the interaction mechanism of ions, 
mediated by the exchange of spin-one particles, such as virtual photons, vector mesons, and gluons. For sufficiently 
large electric charges of the ions, virtual-photon exchanges will eventually play the dominant role. Actually, the effect 
of the replacement a — > Za for the case of charged heavy ions, e.g. for Pb-Pb collisions, exceeds the corresponding 
QCD contribution for typical regions of momentum transfer, where a s ~ 0.1 — 0.2, which seems to be essential in the 
experimental set-up. We shall consider processes of the creation of one and two vector mesons, such as w, J/ip, p, or 
ortho-positronium: 

Y 1 {Z U P 1 ) + Y 2 {Z 2 ,P 2 ) -> V(e,r) + Y 1 (Z 1 ,Pl) + Y 2 (Z 2 ,P^), 

Yi{Z u P x )+Y 2 {Z 2 ,P 2 ) -> V{ ei , ri ) + V(e 2 ,r 2 ) + Y 1 (Z 1 ,P[) + Y 2 {Z 2 ,P 2 -). (1) 

Here, Pj, P! are the four-momenta of the incoming and scattered ions, and e, and r, r t are the polarization four- vectors 
and four-momenta and of the created vector mesons, which obey the transversality conditions e(r)r = ei(ri)ri = 0. To 



, — i 

X 



'Electronic address: ahmadov@theor.jinr.ru 
1 Electronic address: knichl@desy.de 
t Electronic address: scherbak@mail.desy.de 
§ Electronic address: kuraev@theor.jinr.ru 



2 



describe the peripheral kinematics, it is convenient to introduce the light-cone four-momenta pi as linear combinations 
of the incoming- ion four-momenta Pf. 

p 1 =P 1 -\P 2 , P2 = P 2 -vPi, pl=P 2 2 = 0, Pi=m\, P%=ml 2P X P 2 « 2 Pl p 2 = s » m 2 , (2) 

where are the masses of the ions. In the calculation of the differential cross section, the effects of the off-mass- 
shell-ness of the exchanged photons and vector meson must be taken into account. Our approach is based on taking 
the constituent quarks to be QED fermions. The additional factors for the QED amplitudes associated with the color 
and charge of the quarks will be discussed later. 

The details of the wave function of the bound state were discussed in Ref. [2j . For our approach, only one structure 
R is relevant. The virtual-photon polarization effects in the process j*gg — > ip were considered in Ref. [3(. The analysis 
of the inclusive annihilation of heavy quarkonium beyond the Born approximation of QCD was presented in Ref. 
In our work, we shall obtain the differential cross sections for the creation of one or two vector mesons considered as 
bound states of the relevant quarks. 



II. MATRIX ELEMENTS OF THE 2^3 PROCESSES 



To lowest order in perturbation theory, there are two sets of Feynman diagrams involving three virtual-photon 
exchanges (see Fig. 1). The contribution of each of them to the total cross section has the form a ~ ao(aL 2 + bL + c), 
with L = ln(s/M 2 ) being the "large" logarithm. The interference of the relevant amplitudes only contributes terms 
devoid of the "large" logarithm. Below, we restrict ourselves to the consideration of just one of the amplitudes, which 
corresponds to exchanges of one virtual photon with one ion Y\{Z\, P\) and of two virtual photons with the other ion 
12(^2, Pi)- Using the prescriptions proposed in Ref. [l[ for evaluating the matrix element of the peripheral process of 
single vector meson production, we obtain 



M 



' x Y 2 ^Y x Y 2 v _ (47raZi)(47raZ 2 ) 2 / 2 



d 4 q 



(2tt)V3§ 



s 2 N 2 sFC v , 



(3) 



where the factor C v accounts for the color and charge of the quarks and N±, N 2 , and F are given by the following 
expressions: 



1 



JVi = -u(i^)A«(i\), 



1 

2s- 
I MAI 



N2 = T^u(Pi 



F 



s 2 4 



„ Vi - q + m 2 A p 2 - q 3 + m 2 „ 
Pi 7 To 2 Pi + Pi 7 2 Pi 

(P2 - qr - m 2 \P2 - <73r - ^2 

Ti0^ x {p + M 1 )ep^p 2v p 2X . 



u(P 2 ), 



(4) 
(5) 



Here, M is the mass of the created vector meson, and q$ = q 2 — q. The coupling constant A, which we shall 
specify below, measures the strength of the interaction of the vector meson with incoming photons. The quantity 
u{q)0^ vX v{q+) is the matrix element of the subprocess 37 — > qq depicted in Fig. 1. So, we have 



Pl^P 2 uP 2 \ = Pl- 



Qi 



1 . , „ . s . 1 . , . 
-f^-P2(-q+ +Q3 + m)p 2 + —p 2 {-q+ 

L> 2 



-prfciq- ~Q + m)P2 + ~^-p 2 (q- - <73 + m)p 2 

l->2 JJ3 



q 2 + m)p 2 
q+ + qi 



m 



-Pi 



1 



D 2 D, 



Di 

\p2(q- - q + m)pi(-q + + q 3 + m)p 2 + p 2 (q- - §3 + rn)pi(-q + + q 2 + m)p 2 ], (6) 



where q± are the four-momenta of the created quarks and the denominators are given by the expressions D\ = 
(<7_ — <7i ) 2 — m 2 + iO, D 2 = (—q~ + q) 2 — m 2 , and D 3 —[—q + + q 3 ) 2 — m 2 . Following the rules for the construction of 
the matrix element of the quark-antiquark bound state [5[, we must put q + = q_ = p/2 in this expression. Let us now 
introduce the Sudakov parametrization of the loop momenta and the momenta of the quarks and virtual photons: 



qi = P1P1 + qi±, 
q± = a±p 2 + /3±pi 



qz = U3P2 + P3P1 + q3±, q = ap 2 + ppi 



5J 



q±±, 



d q = —dadfid q. 



(7) 



From four-momentum conservation and the on-mass-shell conditions of the quarks, we have 

/?+ = /3- = ^@i, a + a 3 = a 2 ; 

qi+q 2 = qi+q + q 3 =p, sa± = ^—(f + Am 2 ). (8) 

Zpi 

The expressions for the denominators can now be rewritten as: 

Di = -tf+pqi -\(f + 4m 2 ) = + q% + M 2 ) = ~R, 

D 2 = -<f +pq - ^sfoa + iO, 

D 3 = -qi+pqz-\s{a-a 2 )(3i+iQ. (9) 
Performing the integration over the component f3 2 of the loop momentum, we obtain for N 2 



oo 

i J dpN 2 = ^-u(P^p lU (P 2 ) J sdp 



(-2ni)^u(P^u(P 2 ) = -2mN' 2 . (10) 



_sf3 + a + i0 -s/3 + b + iO 

— oo 

Summing the squared moduli of Ni and N' 2 over the spin states, we obtain 

]T|AM 2 = 2, ]T|^| 2 = 2. (11) 
In a similar way, the integral over the Sudakov variable a gives us 

oo 

I' da (-k + -k) = < - 2 " ) h (12) 

— oo 

As a result for the contribution of these terms to the matrix element, we have 

2 (-0^-\2 -i 7i/r A r r ]2 



8(47raZi)(47raZ 2 ) 2 (2Tri) 2 1 Mi A f d 2 q 1 



with the trace 



IT - WlW W5T^T/ ^i Trft (13) 



^TrQ = ^Tv\p 2 (~p + qi + m)pi - pi(^p - qi + m)p 2 ](p + M)e 

= ^Tr[p 2 giPi + Pi qip 2 ]Me = sM(qe). (14) 



The contribution of the last two terms may be found using the relation 



/ 



"'"ok = { - 2 " ] J[3- < 15) 



where D = R — Aqq 3 . The calculation of the relevant trace leads to 

^ Trp 2 {(q-± - q 2 ± + m)(-q+± + q 3 ± - m) + {q~± - q 3 ± + m)(-q+± - q 2 ± - m)](p + M)e 
1 M 

= -2m-Trp 2 gipe = -s— Pi(qie). (16) 

Here, the expression contained within the square brackets is 2m(q 3 + q 2 — q + — g_) = —2mqi, the replacement 
(p + M) — » p is due to the relation pe — 0, and p 2 p — 2 s(3i. Furthermore, we have 

1 1 _ Aq{q 2 - q) 

D R RD ■ 1 ' 
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We thus obtain 



daF = 



(27ri)8gg 3 (egi) MA 
RD 2~~* 



(18) 



It is worth noting that the right-hand side of this equation vanishes in the the limit where the transverse momenta 
of each of the three virtual photons vanish, which is a consequence of gauge invariance. 
The final result for matrix element corresponding to the Feynman diagram in Fig. la is 



M Y 1 Y 2 ^Y 1 Y 2 \ 



M 2 / ff2 y |a3 ^(fle]ij( g 1,fl)C y , 
q 1 ri 



(19) 



with 



J{qi,R) 



d 2 q (qq 3 ) 



where R = q\ + gf + M 2 and q 3 



Q2 - q- 

As shown in Appendix A, the function J has the form 



(20) 



l 



2q 2 - R 



91 R 1 R ~ti 2 
21n ^~ ln -.2 
Z Q 2 92 



1 



qi + M 2 



(21) 



where 



x - 1 



■In 



(1 + ^) 2 
4x 



«1,2 



92,1 



95,2 



M 2 ' 



(22) 



III. DIFFERENTIAL CROSS SECTION OF SINGLE VECTOR MESON PRODUCTION 



The Feynman diagrams contributing to the process of single vector meson production are shown in Fig. 1. The 
phase space volume has the form 



dr, = 



(2tt) 4 d 3 P{ d 3 P{ d 3 P{ 



(2?r) 9 2E[ 2E[ 2E[ 



^5\P 1+ P 2 -P[-P^p). 



(23) 



We introduce the factor d i qi8 A (pi — ch — p[) d 4 q 2 S 4 (p 2 — q 2 — p 2 ) m terms of Sudakov variables and use the standard 
Sudakov parametrization of Eq. ([7]) |l|. Allowing for the vector meson to be unstable, the three-particle phase space 
volume then becomes 



4s (2tt) 9 tt 2 /3i 

where p 2 — sa 2 j3i — (<fi + q?) 2 and we replace the delta function by a Breit-Wigner resonance, 



(24) 



1 



6(p z - M ) -> R{p z ) = - 



TM 



tt (p 2 - M 2 ) 2 + M 2 T 2 ' 



(25) 



where M ss 2m q and T are the mass and the total decay width of the vector meson resonance and m q is the mass of 
the bound quarks. The quantity p 2 = (Pi + P 2 — P[ — P 2 ) 2 may be associated with the missing mass in the process 
of single vector meson production. 

With the matrix element of Eq. ([19]) . we now have 



R 



where 



Z 2 (qle) , . ^i(gle) , . 



<?M+M 2 ) 



qt(qt + M 2 ) 



(26) 



(27) 
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The relevant differential cross section is 



daW 2 6 7r(Z 1 Z 2 ) 2 a 6 A 2 M 4 - d 2 qi d 2 q 2 2 d/3i n . 2 . 

-Z - —R\p ). 



dp 2 



R 2 



Pi 



(28) 



The interference term in Z 2 is canceled by the average over the azimuthal angle. For the case of extremely small 
transverse momentum \q\ \ of ion Y\, a modification of this formula is necessary, which consists in replacing 



1 dpi 1 - Pi dpi 

(qj) 2 Pi ^ {$+mlft)* Pi ' 



(29) 



and a similar replacement for small values of \q 2 \. This leads to the so-called Weizsacker-Wiliams enhancement factor 
in the total cross section: 



^ = a R(p 2 )[Z 2 (L 2 5Li) + Z 2 (L 2 5L 2 ) + c(Z 2 + Z 2 ) + 0(^-, ^)], 



(30) 



where 



32n(ZiZ 2 a 3 ) 2 A 2 



M 2 



(l-ln2), 



ba- 



rn 



1,2 



dx 

x 



5 + x 



2(1 -xf 2 \{l -x) 



l + 2x 



1 In- 



10.4565. 



(31) 



IV. DIFFERENTIAL CROSS SECTION OF TWO VECTOR MESON PRODUCTION 

In the case of two vector mesons in the final state, we must consider the two mechanisms due to the sets of Feynman 
diagrams indicated Figs. 2a and 2b. 

The phase space volume of the four-particle final state is 



, r (27T) 4 3 dp! dp d 2 qid 2 cfd 2 q 2 a a 3 2 1 

4 = (2^ 7r XT 7 ^ dridr 2 R( ri )R(r 2 ) — , [M) 



where r\ = sa'Pi — (qi — (f ) 2 and r 2 — sa 2 (3' — (q 2 + <f) 2 . Also here, we assume that m\ 2 / ' s <C P' <C Pi <C 1. Let us 
introduce another auxiliary four- vector: q' = a'p 2 + P'pi + q' ± . 

The matrix element of the process of two vector meson production mediated by a virtual vector meson (see Fig. 2a) 
has the form 

m y^ymvm = s i^m^l^ (33 ) 

q + My 

where T is given by the expression 

rp = J(qlRl)J(qlR2) + J(q 2 i,Ri)J_{qlR 2 ) 

RiR 2 RiR 2 

with Ri = q 2 + q\ + M\ , R 2 = q 2 + q\ + Mf , R x = q 2 + q\ + Mf, and R 2 = q 2 + q\ + Mf . Here, M v is the mass of 
the virtual vector meson, g is its coupling constant, and Mi t2 are the masses of the created vector mesons. 
The cross section of the process of two vector meson production is 

da { a ] 2 6 a s (MiM 2 ZiZ 2 ) i 4 \2 r>/ 2\ n/ 2\\rp\ 2 / -*-* \ 2 / \2 //~fVi s~iV2 \2 d 2 qi d 2 q 2 d 2 q dPi dj3 
- ~ — - g^{AiA 2 ) i R{r{)R{r 2 ^\T\ z {qeiY{qe2Y{C 2 1 C 2 2 Y — -. (35) 



dr 2 dr\ ~° {(P + M 2 ) 2 y " l ' u " V2 ' M ^ L > ^ ^ 2 ~ 2 > tt3 Pi p- 

After integration over d 2 q, we obtain: 

do^ _ 8ZlZjgW .dPidP 2 ( M V M V \ 2 v v d^ 2 qj f 2 

-jprfiS- (l + 2cos 6i 2 ) — -{AiA 2 ) [-^^) R(ri)R(r 2 )(C 2 C 2 ) -^P^,^), (36) 



G 



where cos#i 2 = coseie 2 and 



M^Ml> J (q^ + M 2 ) 2 Ml 
o 

is evaluated in Appendix B yielding the numerical values presented in Table I. 

For the case of two-gluon exchange (see Fig. 2b) , the matrix element has the form 

rYiY^YiYMV* _ ■ 2 12 tt 3 (Z 1 Z 2 a 2 a 2 ) (M Vl My 3 ) 2 N,N 2 A A , 



M t r '^ r ' f ^' Va = is v 1 ' % y ' V1 V2/ — A 1 A 2 QC^Cp, (38) 



where 



I(Ri,R 2 ) ^ . I(Ri,R2) -w_ -,s , Qn s 

<2 = p p — (eigi)(e 2 £ftj) + p p — (eig 2 )(e 2 gi). (39) 

The quantity I{R\, R 2 ) is given in the Appendix A. For the cross section, we may write 

d °l 2) __ ^MZiZ 2 a 2 air(M Vl M V2 ^ 2 d 2 q 1( Pq 2 d?q . v v - a 
*prf (fl? RirMrJQ — — (AiA 2 ) ^ (C 2 C 2 ) . (40) 

Performing the integration over d 2 cf, we obtain 



$ + $ + 2cos 2 6 12 G 



W, (41) 



where the functions $, and G are evaluated in Appendix B yielding the numerical values presented in Table II, 



W= d 4 d 4 (42) 



is the Weizsacker- Williams enhancement factor, and 



Xl - 2 = m^mV 2 - (43) 

,™ m „of ^1„™ ^2 > ^2 i «2^ 2 x<2 , =2 , „, 2 ^ 2 



In the region of small values of qf 2 , we must replace qf — > qf + /3 2 m 2 , q 2 — !• <z 2 + a 2 m 2 

If, 2 + i«„ ; - | 



We note that the interference term of the two amplitudes is absent, so that \Mf' + M^\ 2 = |M b ^| 2 + \M^\ 2 . 



V. CONCLUSION 

We studied the production of one or two vector mesons in peripheral heavy-ion collisions at high energies. In the 
case of Za > a s , a simplified version of a general theory Q can be used to lowest order of QED and QCD that is 
based on the subprocesses 7*7*3 — ► V and "f*gg — > V. In our study of two vector meson production, a vector meson 
can also appear virtually as an intermediate state. In this case, it is important to replace it by the relevant vector 
reggeon state, with the Regge trajectory 

a v ((q) 2 ) = a v (0) + a' v (0)(q) 2 « a v (0) « 1/2. (44) 
This results in an additional factor 

M ~^§, (45) 

so J P 

where sq ~ 1 GeV and p 2 is the missing mass square, in the cross section da^ . 

When constructing the invariant mass square of the decay products of one vector meson, p 2 = (Pi +P 2 — P[—P 2 ) 2 = 
sf3ia 2 — (qi+q 2 ) 2 , we take into account that the part sj3\a. 2 is the combination Ei) 2 — (^2pi Z ) 2 , with Ei andj?i 2 being 
the energies and z components of the three-momenta of the decay products, while the part — (<fi+<f 2 ) 2 is the contribution 
from the transverse components — {^2Pi±} 2 - Here, it is understood that the z direction is taken along the beam axis 
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in the center-of-mass frame. Such is the case for two jet production with r\ = (Pi — P[ — q 1 ) 2 = sa'ftx — (<fi — q 1 ) 2 
and r 2 = (P 2 -P>+q') 2 = S a 2 f3'-(q 2 +<f) 2 . 

The coupling constant A of the meson-photon interaction in the case of single vector meson production, appearing 
in Eq. (5), is given by A — a 3 ^ 2 /(2^/tt) for ortho-positronium and by Ai = 2fyJMy i , with f p = f u = 0.21 GeV and 
f^p = 0.38 GeV, for the w, p, and J/tp mesons, respectively (sec Rcf. [6] for details). 

In the case of single vector meson production, the color and charge factors are C v = 3 ^ Q 3 = | for the p and ui 

u, d 

mesons, and C v = | for the J/tp meson. In the case of two vector meson production through the mechanism shown 

in Fig. 2a, we have Cvg = 3 Q\ = § f° r the p and co mesons, and C 2 X = | for the Jjip meson. In the case of the 

mechanism shown in Fig. 2b, we have C^ 1 — Sj^Qq = 1 f° r the p and lo mesons, and C^ 1 = 3| = 2 for the J/ip 
meson. 

Note that the mechanism involving single 7* exchange (see Fig. 2b) yields a L 4 enhancement, whereas double 7* 
exchange (see Fig. 2a) only produces a L 2 enhancement. For pp collisions at the LHC, the "large" logarithm is as 
large as L = In -^7 ~ 7. 

We do not consider gluon exchange between heavy ions and vector mesons to avoid channels with ion exitation. 



Appendix A 



In this section, we shall explain how to calculate the two-dimensional Euclidean integrals appearing in Eqs. (fl9|) 
and (1391). 



d 2 k k(q 2 - k) 
71" k 2 (q 2 - k) 2 D'' 



I(R 1 ,R 2 ) = 4 



d 2 k (k{q-k)f 
v k 2 (q-k) 2 D 1 D 2 



(Al) 



where D = R- 4k(q- k), R = qf + q\ + M 2 , A = Ri - 4fe(g - k), and R { = q 2 + q 2 + M 2 . For the first one, we have 



4J 



with 



Ji = 



d 2 k 



1 



71 (0?2 - k) 2 + X 2 )(k + \ 2 )D 
Using Feynman's trick of joining the denominators, 



d 2 k R-D 

— =; = hm 

K k 2 {q 2 -k) 2 D A^O 



Jo = 



RJi — Jq 



d 2 k 



(k 2 + \ 2 ){{q 2 - k) 2 + A 2 ) 



(A2) 



(A3) 



— = I dx 



ab 



(xa + xb) 



2 ' 



(A4) 



with x = 1 — x, we obtain for Jq: 



Jo = dx 



d 2 k 



* [(k - qVx) 2 + qVxx + A 2 ] 2 



dx 



q\xx 



A 2 gf A 2 ' 



(A5) 



For Ji, we have 



1 i 
4Ji = dx 



ydy 



[q^x^ + ^R-y^+yX 2 } 2 ' 

where y = 1 — y. Introducing the variable t = 1 — 2x, we may cast this into the form: 

1 1 
Ji = 4 / ydy 



dt 







(A — Bt 2 ) 2 ' 



(A6) 



(A7) 
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where 

A = f 2 y 2 + y(R-qi 2 y)+Ay\ 2 , B = cf 2 y 2 . 
Performing the integration over t, we obtain 



(A8) 



Ji 



ydy 



a 3 / 2 b 1 / 2 



2(ABW 2 a 1 / 2 +b 1 / 2 

1 In 



A - B A 1 / 2 - B 1 / 2 



(/ 1+ / 2 )Q 2 , 



(A9) 



where 



h=2 



ydy 



[vip-V) +4yo-]T' 



dy , T 1 ' 2 + y 
In ■ 



y3/2 JU/2 



(A10) 



with 



T = y 2 +y(p-y), P=^>' CT 



'7 2 



(All) 



The first integral I\ contains an infrared singularity. Introducing the small parameters a and e, with u < 6 < 1, 
we rewrite it as 



l-e 



h = 2 



/ 



ydy 



+ 2 



y(p-y)[p-i + y{2-p)} 

l-e 



dy 



yp + 4a 



— In - 

P 4ct p{p-l) 



[(p-l)ln(p-l) + lnp]. 



(A12) 



For the second integral I 2 , which is infrared finite, the substitutions T = t z , y = ■ L -j^, a = p — 1, and a + b = 1 yield 



/a = 



2 } dt tb + t 2 - a _ 2 f dt(l - t) 
b J T 2 U tb- (t 2 -a)~ b 

1 



b + 2t 



b-2t 



tb + t 2 - a tb-t 2 +a 



4 / dt- 



t z + a 



t(l + t)(t 2 -a 2 ) b 



2 [lna - 2 In 2 + ln(a + 1)]. 



The total answer for Ji is 

Ji(p) = 2 

For the sum J = RJ\ — Jo, we obtain 



V 1 



P cr p(2- p) 



1 , i? 2 

In 



(21np- 2 In 2 — ln(p — 1)) 



1 , (x + 1) 2 
In ■ 



2q 2 -R 4(R-c%)(% gf + M 2 x — 1 4a; 



where x = gf /(gf + M 2 ). 
For the integral /, we have 



4/ 



d 2 fc 1 (R 1 -D 1 )(R 2 -D 2 



Jo 



R 2 



-Ji(Ri), 



7T k 2 (q-k) 2 DiD 2 R1-R2 R1-R2 

where Ji given in Eq. (|A9[) and i?i j 2 = g?2+<7 2 + -^12- The infrared singularity is canceled using 

1 



i?2 , R 
■ In 



Ri , -R 2 

In ■ 



i?i - i? 2 pq 2 - R 2 4q 2 (R 2 ~f) 2<f - Jfr 4<f(i?i-g 2 ) 



(A13) 



(A14) 



(A15) 



(A16) 



(A17) 
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For the case of Ri = R 2 = R = q\ + q + M 2 , we have 



d f d 2 k k(q — k)(R - D) 



7T k 2 (q-k) 2 D 



A—RJ(R,<f)=A 



2<Z? 



In- 



R 2 



R-q{ 2qJ-R Aq({R 2 -q{) 



(A18) 



For the cross section of single vector meson production, integrated over q 2 and q 2 , we have 



d P 2 R{p 2 ) 



Using 



\^(Z Y Z 2 c? ) 2 A 2 
M 2 



dx 



■In" 



(x + If 
Ax 



(X + 1) 2 (X- l) 2 

u 

to? m 2 , . 1 

(A -► a 2 , Z x o Z 2 , Pl = j± -> P 2 = 



dx 



(x 2 - l) 2 

we obtain the expression given in Eq. (31). 



ln 2 ( 



(■^ + 1) S 
Ax 



dt ■ t 



2(1 -ln2), 



(A19) 



(A20) 



Appendix B 



In this section, we shall explain how to evaluate the integrals appearing in Eqs. (|36[) and (|41j) relevant for the 
mechanisms based on vector meson (see Fig. 2a) and two-gluon (see Fig. 2b) exchange, respectively. 
In the first case, we have 



P(xi,x 2 ;pi,p 2 ) 



dx • x 



(X + l) 2 ' 



(Bl) 



with 



where 



i(x 1 ,r 1 )i{x 2l r 2 ) i(xi,ri)i(x 2 ,r 2 ) 



r\r 2 



r\r 2 



i(x, r) 



■In 



2.T — r Ax(r — x) ' 



(B2) 



(B3) 



and n = x + x\ + pi, r 2 = x + x 2 + p 2 , r x = x + X\ + p 2 , r 2 = x + x 2 + p x , pi = M Vl /M v , and p 2 = M 2 2 /M v . 
Numerical values of P(xi,x 2 ; 1, 1), appropriate for the important case pi = p 2 = 1, are listed in Table [IJ 
On the other hand, we have 



$(xi,x 2 ;pi,p 2 ) 
®(xi,x 2 ;pi,p 2 ) 
G{xi,x 2 ;pi,p 2 ) 



dx 



dx 



nr 2 

j{ri,f 2 ) 
fif 2 



= I dx 



J( ri,r 2 )j(ri,r 2 ) 
r\r 2 fi,f 2 



(B4) 



where 



j(ri,r 2 ) = 



n - r 2 



r 2 



■In 



ri 



In- 



2x — r 2 Ax(r 2 — x) 2x — ri Ax(ri — x) 



(B5) 
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U.UUUU40 


A AAAAQQ 

U.UUUUoO 
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U.UU4 1 
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U.UUUU4y 


A AAAA1 
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Q act in-b 
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O JO 1 A t> 

z.4o • 1U 
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1 

1 


A AAA 1 A 
U.UUU49 
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U.UUUU49 


A AA 1 A 1 A~ti 


Z.o4 ■ 1U 


7 or i A — 7 
(.00 • 1U 


o p;o i a™ 7 
z.oz ■ 1U 


O 1 Qfi 1 A — "( 
Z.loO ■ 1U 


O AQ 1 A - 1 

z.Uo ■ 1U 


1.0 


A AAAAQO 

U.UUUUoz 


A AAAA1 

U.UUUU1 


O Q A 1 A - 6 

z.o4 • 1U 


1 OA 1 A - 

l.o\) ■ 1U 


O 7Q 1 A - ^ 

z. • 1U 


AR 1 A - S 

0.40 ■ 1U 


Q Q K'V 1 A - S 

o.oo/ ■ 1U 


Z.00 ■ 1U 


2 


0.00004 


3.66 ■ lO -6 


7.85 ■ 10~ v 


2.78 ■ 10~ Y 


1.29 ■ W 


4.59 ■ 10"* 


2.52 ■ 10"* 


1.76 ■ 10" 8 


3 


0.000046 


2.51 • 10" b 


2.52 ■ 10" Y 


6.46 • 10" 8 


4.59 • 10~ 8 


4.02 ■ 10"* 


3.21 ■ 10"* 


2.45 • 10" 8 


4 


0.000045 


2.48 • 10~ b 


2.14 ■ 10~ Y 


3.36 • 10~ 8 


2.52 • 10" 8 


3.21 ■ 10"* 


2.91 ■ 10"* 


2.35 ■ 10~ 8 


5 


0.000038 


2.23 • 10~ b 


2.029 ■ 10"' 


2.85 ■ 10~ 8 


1.76 ■ 10 -8 


2.45 ■ 10~* 


2.35 ■ 10"* 


1.96 • 10" 8 



Table I: Values of the function P(x\, x%\ 1, 1), denned in Eq. (|B1|) . for different values of x\ and Xi- 



XI \ X-2 


0.1 


0.5 


1 


1.5 


2 


3 


4 


5 


0.1 


1.67 


1.006 


0.623 


0.428 


0.314 


0.191 


0.129 


0.094 


0.5 


1.006 


0.6142 


0.386 


0.268 


0.198 


0.123 


0.08399 


0.0615 


1 


0.623 


0.386 


0.246 


0.173 


0.129 


0.0808 


0.0559 


0.0414 


1.5 


0.428 


0.268 


0.173 


0.1221 


0.0918 


0.0583 


0.04078 


0.0304 


2 


0.314 


0.198 


0.129 


0.0918 


0.0695 


0.0446 


0.0314 


0.0235 


3 


0.191 


0.123 


0.0808 


0.0583 


0.0446 


0.029 


0.02073 


0.0157 


4 


0.1297 


0.08399 


0.0559 


0.04078 


0.0314 


0.02073 


0.149 


0.01142 


5 


0.094 


0.0615 


0.0414 


0.0304 


0.0235 


0.0157 


0.01142 


0.0088 



Table II: Values of the function &(xi,X2\ 1, 1) = &(xi,X2\ 1, 1) = G(xi,X2; 1, 1), defined in Eq. (|B4Jl . for different values of xi 
and X2- 



In the case r\ = r%, we have (see Eq. (|A18|I ) 

8x ,/(x + l) 2 \ q? . , 

j(n,n) = 4 - j—^ Xn y—^j > * = 

For pi — p2, we have <fr(xi, xi\ pi, p{) = X2', pi, pi) — G(xi, X2', pi, pi)- Numerical values for the important choice 

pi = p2 = 1 are listed in Table [TT1 
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Figure 2: Feynman diagrams pertinent to two vector meson production in peripheral heavy-ion collisions via intermediate a) 
vector meson and b) two-gluon states. 



